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We study the creation of nonlocal states with ultracold atoms trapped in an optical lattice. We show that these 
states violate Bell inequality by measuring one- and two-body correlations. Our scheme only requires beam 
splitting operations and global phase shifts, and can be realized within the current technology, employing single¬ 
site addressing. This proposal paves the way to study multipartite nonlocality and entanglement in ultracold 
atomic systems. 
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Bell’s theorem Cl unveils the incompatibility between 
quantum mechanics and a generic local hidden variable the¬ 
ory via a set of inequalities that can be experimentally tested. 
The notion of Bell nonlocality (see (2) for a recent review) has 
been originally proposed and mainly developed in the case of 
two parties. It is now recognized as a crucial resource for 
several quantum information processing, as secure commu¬ 
nication El IS and certified random number generation 0. 
Violations of Bell’s inequalities have been reported in a va¬ 
riety of bipartite physical systems, giving a strong evidence 
that nature is nonlocal 0. Experiments have been mostly 
performed with photons [7- 1T| and, in recent years, also with 
ions MM and atom-photon hybrid systems llT5lfT6ll . Quan¬ 
tum nonlocality in many-partite systems is, by far, more com¬ 
plex and less developed than in the bipartite case. One crucial 
obstacle is that multipartite Bell’s inequalities [ 171 HI usu¬ 
ally involve high-order correlations that should be measured 
for N parties in order to demonstrate nonlocality. This poses 
serious experimental - but also theoretical - challenges for 
large N. Presently, the experimental violation of multipartite 
Bell’s inequalities has been achieved with three lfl91-f2D and 
four (22| [23) photons and trapped ions (24), whereas most 
of these are limited to GHZ states. Much less is known about 
the possibility to investigate quantum nonlocality for ultracold 
neutral atoms. These systems offer the practical advantage of 
large detection efficiencies and a variety of control techniques, 
including single site addressing and fine manipulation of in¬ 
ternal and external degrees of freedom. Recent progresses 
have evidenced that ultracold atoms are optimal candidates for 
the creation of entangled states with a large number of parti¬ 
cles 125 -291. However, while entanglement is generally a pre¬ 
requisite for the violation of Bell’s inequalities (30) . to date, 
there are only few proposals (3X1433) to demonstrate quantum 
nonlocality in these many-body systems. 

Here we propose a realistic experimental protocol to create 
and observe Bell’s nonlocality of an arbitrary large number 
of neutral atoms trapped in an homogeneous one-dimensional 
optical lattice. The crucial ingredient of the proposal is the 
coupling between hyperfine states of atoms sit at neighboring 
lattice sites. This coupling can be accomplished by resonant 
light interaction and a moving state-dependent lattice poten¬ 
tial, as outlined in Fig. [I] We notice that the combination of 
internal Rabi splitting and the coherent transport of atoms in 



(a) Mott state preparation 



(b) nearest neighbor coupling 





(d) on site coupling 


FIG. 1. (color online) Quantum nonlocality of neutral atoms in an 
optical lattice (here limited, for illustration sake, to few lattice sites). 
The protocol consists of four operations involving internal (hyper¬ 
fine) states and it is implemented by a state-selective optical lat¬ 
tice. (a) Preparation of a system of two-level atoms (a and b indi¬ 
cate hyperfine levels, l the lattice site) with one atom per lattice site, 
each atom being schematically represented by a dot. (b) Balance 
pulse coupling neighboring wells: this is made of a 7t/2 Rabi pulse, 
schematically represented by a green waving line, (b.I) followed by 
a relative shift of the state-selective lattice (b.II). (c) Collective phase 
shift, 0 of (p, see text; (d) Balance pulse coupling the internal lev¬ 
els of each atom. The protocol ends with the measurement of the 
number of particles in the internal ground level. 


a spin-dependent optical lattice has been first experimentally 
demonstrated in (34) . This capability has been further used 
to generate entangling gates for neutral atoms (351 36), super¬ 
exchange interaction ED, quantum random walks (38) . and 
spin-squeezed states E7I . Our protocol for the violation of 
Bell inequalities consists of four steps. A Bose gas is ini¬ 
tially prepared with one atom per lattice site, see Fig. [TJa), 
each atom having two internal levels indicated as a and b. 
This can be accomplished by starting from a superfluid gas 
and then raising the potential barriers in such a way to enter 
the Mott-insulator phase (39) . eventually suppressing tunnel¬ 
ing between neighboring wells. The system is thus prepared 
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in the state l^o) = (^) | a) z , where | k) l indicates an atoms in 

i 

the internal state k = a,b at lattice site l = 1, 2 ,..., N. As 
second step, each atom is placed in the coherent superposi¬ 
tion | a) l (| a) l + |6)^ + 1 )/v / 2, see Fig. lib). This operation 
(linear in the N particles) creates correlations between lat¬ 
tice sites, which are our parties. It will be shown that these 
correlations are nonlocal as they are responsible for the viola¬ 
tion of Bell’s inequalities. Nearest neighbor coupling can be 
experimentally realized following two stages: I) a 7r/2 Rabi 
pulse for atoms at each lattice site, | a) l (| a) l + \b)j)/y /2 
[Fig. |TJb,I)] ; followed by II) a shift of the state-dependent lat¬ 
tice potential \b) l | b ) l+1 [Fig. [Tib, II)] that can be realized 
by changing the polarization angle of two linearly-polarized 
and counter-propagating laser fields mm. As third step, we 
apply a local phase shift e l0in9 ’ 1 at each lattice site [Fig.yjc)], 
where n g j is the number operator and 6 i are a set of local 
phases. Below we show that the optimal choice for the phase 
shifts involves imprinting the same phase to each atom. This 
is an important simplification for the experimental implemen¬ 
tation of the protocol. Phase imprinting should happen on 
a time scale much faster than the interaction time scale be¬ 
tween the | g) l and |e) i+1 atoms. We point out that, differently 
from the creation of entangling gates [41, 42), our protocol 
does not involve interaction between atoms. The final oper¬ 
ation is an on-site coupling pulse W)i (I °)i + \b)i)/V 2 
and | b) l (\b) l — | a) t )/V2 [Fig. [Tjd)], formally described 

as e~ l ^ J y with J y = \Y^i=i° l y an( 3 the Pauli matrix for 
the Ith site. The protocol ends by counting the total number 

N 

of particles in the internal ground state N g = ^ n g j, from 

i=i 

which we obtain a dichotomic measurement given by the par¬ 
ity P g = (— l) N f> . We emphasize that the whole protocol does 
not involve vibrational states of the well. 

Our criterion for witnessing quantum nonlocality of the ul¬ 
tracold atomic system is based on the violation of a novel class 
of multipartite Bell’s inequalities, introduced in (43) and in¬ 
volving only one- and two-body correlations. These inequali¬ 
ties are 

B N {9 u <pi) = a(s 0 + |r) + |5oo + y5oi-^ > -47 

( 1 ) 

where 


N 

= e<M°> 

i =1 


N 

Ski = E (2) 


i,3= 1 


denote the one- and two-body correlations, with (k = 
0 , 1 ) representing the two different local measurement observ¬ 
ables realized in well l : 


= g _ l) Ng e * 2 

— g — i(Pl n g,l Jy ( 


(3) 


ferent detector settings in the Clauser-Horne-Shimony-Holt 
(CHSH) scheme (44|. We choose a = N(N — 1) (ff ] - f) 
and 7 = n(n^-i) , * n accor( j ance ^ parameters used in 
l43l to trace nonlocality of Dicke states, for which the clas¬ 
sical limit is set by (3^ = N ^ N 2 ~^ To illustrate 

our proposal, we will first focus on the simple case of two 
atoms trapped in two wells and then generalize the discussion 
to many atoms and many wells. 

Two-wells case. For two and three wells our proposal re¬ 
duces to the scheme first proposed by Yurke and Stoler [45, 
|46l for optical systems. In the double-well case, the final state 
is IV’fin) = (e^ 6 ' 1 + 6 > 2 )/ 2 |' 0 i) + V^ 2 »/ 4 , where 

|Vt) = c (|1010) + |0101>) + is (|1001) + |0110>), (4) 

\^ 2 ) = e i01 (|0200) - 12000)) + (|0002) - |0020>), (5) 

c = cos 6>1 + 6>2 , s = sin 0l +° 2 and \n g ^n e ^n g ^n e ^) repre¬ 
sents the state of occupation numbers in each mode of the sys¬ 
tem, 0 < n gi i < 2. From Eq. 0 , we can obtain analytically 
the probability of the possible measurement results. Similarly 
to Ref. [45], let us denote with P{e \, 62 ) the probability that an 
event e m occurs in the well m. As we consider only bosonic 
atoms, these events are elements of the set [0, G, E , G 2 , E 2 ], 
where G represents one particle in its ground state and E one 
particle in its excited state. We can distinguish three subsets of 
events obtained by the measurement of the population of each 
atomic state. These subset are: firstly the case where both the 
two wells are populated with atoms in the same internal state 
A = {GG, EE}, the second one represents the case where 
both the two wells are populated but with atoms in different 
internal states B = { GE , EG} and finally, we have the case 
where only one well is populated with the two atoms in the 
same state G = {G 2 0, E 2 0, OG 2 , OE 2 }. The corresponding 
probabilities are: 


-P(e i,e 2 ) = < 


jCOS 2 (7i7) 

i sin 2 (£i±^) 

1 
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if € 1 , 6 2 eA, 
if 6 1 ,€ 2 eB, ( 6 ) 

if ci, e 2 £ G. 


We can see here that independently of the local phases, if one 
measures the atomic population in one well, one can predict 
whether the total results belongs to the set A (J B or G. This 
is coherent with local realism as it is a simple consequence of 
the conservation of the total number of particles in the system. 
As a consequence, we can directly ignore the part of the state 
belonging to the set G and test the nonlocality for the reduced 
state |Vh )/2 which we propose to study using the criterion 
0 . We stress here that for two particles the class of Bell’s 
inequalities ([!]) reduces to the well-known CHSH inequality 
(44) : 


—2 < 500 2 — + S 01 < 2. (7) 


The two measurements differ by local phase shifts, 61 and ipi, Using the set of measurement described previously leads to 
applied to the ground state. These correspond to the two dif- the following expression for the mean value of the Bell’s op- 
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erator £2(01,02,^1, ¥>2): 

B 2 (d 1 ,0 2 ,lfii,lfi 2 ) = cos (6 1 + 0 2 ) -cos(<^i + <p 2 ) 

+ cos (ipi + 0 2 ) + cos (Oi + ip 2 ). 

The optimal choice for the phase shifts can be rewritten as 
6 \ + O 2 — (pi + 62 ) — (0\ + ip 2 ) = w and — (pi + 

P2) = 3cj as in Ref. ED- The expression 0 becomes — 2 < 
3 cos(cc) — cos(3cc) < 2 which is maximal for uj = 7r/4 and 
we find B^ rnax>) ~ 2y/2. This violation corresponds to the 
maximal violation predicted by the CireTson’s bound El. it 
is worth to note here that a violation of the CHSH inequality 
is also observed without the reduction of the state for the 
same set of values of 61 , 62 , pi and P2 for which we then find 
a smaller violation i.e. nax ^ ~ 2.41 in accordance with 

Refs. OH and 145). 

Many-wells case. We now generalize the discussion to 
the many-well case. In particular, we emphasize that the 
case N > 2 requires postselection of the output results and 
therefore single-site imaging of the optical lattice, a tech¬ 
nique which has been successfully demonstrated experimen¬ 
tally [48|. In the two-wells case, post-selection of the state 
allows a stronger violation of the Bell’s inequalities but is not 
compulsory to observe it and has not been considered in the 
previous section. After the last beam-splitting operation of 
the protocol outlined in Fig. [T] we obtain (see for instance 
Eq.0 for the double-well case) all wells populated by only 
one atom, with a probability l/2 iV_1 , or at least one empty 
well, with a probability 1 — 1/2 N ~ 1 . In the latter case, ob¬ 
serving the number of particle in N — 1 wells allows to know 
with certainty the population of the last one, independently of 
the local detectors settings {6k, Pk }• As a consequence, the 
correlations involving states with at least one well empty are 
classical correlations. A violation of Eq. 0 requires to ignore 
this case. It should be noticed that the same post selection has 
been considered and proved necessary elsewhere mm in 
order to observe the GHZ contradiction with three particles 
l49l . After post-selection, we can express the mean value of 
the Bell’s operator in Eq. 0 as 

N 

B N (6 k , pk) = E a cos 6k + P cos pk~\~ 
k= 1 
N 

+ 2 cos (^ k + Qi) + S cos (o k + <pi) + 2 cos (v>k + <pi)- 

k,l =1 Z 

kytl 

(9) 

We further minimize this expression as a function of 6k, Pk 
using a genetic algorithm [50). In Fig.[2j we plot the maxi¬ 
mum quantum violation of Eq. ([!]) normalized by the classical 
bound, £jv = / mm.g ktV>k [B N {0 k , ifk)\, as a function of 

the number of parties N of the system. The figure compares 
(crosses) with ^ ob = (3^ / mm e ^ [B N {6,p)\ (orange 
line) obtained with a global optimization assuming that all the 
local phases are equal. We thus conclude that the case 6k = 6 



FIG. 2. (colors online) Maximal relative violation of Eq. ([I), 
optimized over different local phases, as a function of the number of 
parties N (crosses). The orange line is where all local phases 
are equal, 6 k — 6 and p k = p. The dashed grey line is a power law 
fit 1.5/A. 


and pk = p is very close to being an optimal choice of lo¬ 
cal phase. This is of particular interest from the experimental 
point of view as it is far easier to shift the whole system with 
the same phase for each well instead of applying a different 
phase shift to each party. Furthermore, a simple fit of our 
results shows that Bn ~ 1/N, comparable to the result ob¬ 
tained in Ref. |43) for Dicke states. Figure[3]shows the behav¬ 
ior of our relative quantum violation as a function of the phase 
shifts 6 and p for different number of parties involved. There, 
we plot only the values of < 1, i.e. the colored regions 
represent the region where we observe a violation. We see 
here that as the number of parties increases, the robustness of 
the violation against the fluctuations of the phase shifts 6 and 
p increases as well. 

Discussion. Having in sight the possibility to realize exper¬ 
imentally the proposed scheme, some technical issues must 
be studied further. First, let us emphasize that the above dis¬ 
cussion, focused on Bose gases, can be easily generalized to 
fermions. Furthermore, we do not need to assume the parti¬ 
cles originate from a common source (an initial superfluid, for 
instance): interestingly the protocol works even if the particle 
have never seen each others EL0305). 

A possible obstacle to the observation of nonlocality is 
given by interactions between atoms which can affect the first 
beam splitting operations. To evaluate how the interactions 
impact the violation of the inequality 0, we generalize the 
scheme of Fig. [T] replacing the linear beam-splitter operation 
e ~ l 2 J y in Eq. ^ by a non-linear one U = e~ l % (A+xA) 

[5JJ, where J z = \ J2iLi ^ • The coefficient x represents 
the relative strength of interaction with respect to the coupling 
between wells. We then estimate the role played by the inter¬ 
actions on the Bell operator B 2 (6i, 62, Pi,P2) (f° r simplic¬ 
ity we consider only the two-wells case here), for optimized 
phase shifts. This calculation can be performed analytically 
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FIG. 3. (color online) Relative quantum violation oh taking the same local phase shift, 0 and p , at each lattice site. Different panels refer 
to N = 2,4,10,30 (from left to right). The color scale is cut at ^ ob = 1 and violation of Bell’s inequality (£^ ob < 1) is obtained in the 
colored regions. 


for two wells flOll . In the Rabi regime, x ^ 1, we obtain 


B 2 



( 10 ) 


where the first order in x cancels exactly (40l . The maximal 
violation of Bell’s inequalities is decreased by interaction, to 
the second order in x- It is possible to demonstrate [40] that 
this result holds also for more than two wells. This suggests 
that interactions among particles should not be a key problem 
for the violation of Bell’s inequalities in our system. 

Conclusions. Neutral atoms are currently playing a key role 
for the creation many-particles entanglement [ 25-291. Yet, the 
investigation of nonlocality in these systems has been scarcely 
considered. In this manuscript we have proposed a scheme to 
create quantum nonlocality of an ultra cold gas trapped in an 
optical lattice, requiring only linear operations (balanced cou¬ 
pling between internal state of neighboring lattice sites), col¬ 
lective phase shifts, and involving an arbitrarily large number 
of neutral atoms. Nonlocality can be demonstrated via the vi¬ 
olation of a novel set of Bell’s inequalities recently proposed 
in ns and involving only one- and two-body correlations. 
For two wells, the protocol reduces to the one first proposed 
by Yurke and Stoler in E3. For many wells, the relative vi¬ 
olation goes as 1/N as a function of the number of atoms in¬ 
volved in the measurement. The violation of the Bell inequal¬ 
ity requires postselection, as typical [ 31 ,46]. Our work paves 
the way toward the realization and demonstration of quantum 
nonlocality with ultracold atoms - involving a large number 
of parties - and their use for quantum information protocols. 
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